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Abstract. We analyze the asymptotic behavior of the Apostol-BernouUi 
polynomials Bn{x; A) in detail. The starting point is their Fourier series 
on [0, 1] which, it is shown, remains valid as an asymptotic expansion 
over compact subsets of the complex plane. This is used to determine 
explicit estimates on the constants in the approximation, and also to 
analyze oscillatory phenomena which arise in certain cases. 

These results are transferred to the Apostol-Euler polynomials £n{x; A) 
via a simple relation linking them to the Apostol-BernouUi polynomials. 

1. Introduction 

The family of Apostol-Bernoulli polynomials Bn{x; A) in the variable x 
and parameter A € C were introduced by Apostol in [1], where they are 
defined by means of the power series expansion at of the meromorphic 
generating function 



(1) 9ix,X,z) = j^^—j = ^Bn{x;X) 

n=0 

The set S of poles of g is {lirin — log A : n € Z} when A ^ 1 and {27rin : 
n € Z, n 7^ 0} when A = 1, due to the fact that is a removable singularity 
in the latter case. Here and throughout this paper, log is the principal 
branch of the logarithm; namely, for A 7^ 0, log A = log|A| + iargA, with 
— vr < argA < vr; in particular log 1 = 0. This change in the set of poles 
is reflected in various discontinuities as A — ?• 1; for example, the radius of 
convergence of the series in (1) is 27r when A = 1 and [ log A| when \ ^ 1. 

The value A = 1 corresponds to the classical Bernoulli polynomials, i.e. 
Bn{x; 1) = Bn{x), but it is certainly not the case that Bn{x) = limA_!.i Bn{x; A). 
There is a limiting relationship between 0„(x;A) and Bn{x) as A — > 1, but 
it is not immediately obvious. Another aspect of this discontinuity is that, 
although Bn{x) is monic of degree n, for A 7^ 1 the degree of Bn{x; X) is 
n — 1 and its leading term is n/(A — 1). 

We follow the convention of denoting Bn{X) = S„(0;A); these are the 
Apostol-Bernoulli numbers. In fact, Bn{X) is a rational function in A with 
denominator (A — l)""*"^ and coefficients related to the Stirling numbers. The 
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classical Bernoulli numbers are given by i?„ = ;B„(1). The above remarks 
concerning the discontinuity at A = 1 also apply here. 

The case A = is trivial; indeed Bo{x; 0) = and 0) = —nx"~^ for 

n > 1. In particular Bi{0) = —1 and Bn{0) = for n ^ 1. For this reason 
we will assume A 7^ in what follows. 

Dilcher showed in [3], using properties of the Riemann zeta function, that 
the Bernoulli polynomials satisfy 

hm ^ '.^\,^ B2n{z) = cos(27rz), 
n^oo 2(2n)! 

1™ L \ iM B2n+i{z) = sin(27rz), 

uniformly on compact subsets of C. In addition, the difference between the 
nth term and its limit is found to be of the order 0(2~"), with the implicit 
constant depending exponentially on \z\. The authors showed in [5] how 
these facts also follow easily, at least on [0, 1], from the Fourier expansion of 
the Bernoulli polynomials, including the quantitative bounds for the differ- 
ences, and also precise estimates on the rate of convergence, namely bounds 
for the successive quotients of these differences. 

The purpose of this article is to obtain analogous asymptotic estimates 
for Bn{z; A), valid for any z € C. In short, the central result of this paper 
is that the Fourier series (4) of Bn{x]\) for x € [0, 1], which a priori repre- 
sents it only on this interval, is actually valid on the entire complex plane as 
an asymptotic series representing Bn{z;X) for z G C. From this we deduce 
explicit asymptotic estimates for the Apostol-Bernoulli polynomials which 
include the pattern mentioned above, namely, geometric order of decrease 
for the differences between Bn{z; A) and its asymptotic approximations, with 
implicit constants that are exponential in as well as estimates for suc- 
cesive quotients of these differences. The behavior of the approximations 
varies considerably depending on A, with A = 1, studied in [5], turning out 
to be the exception rather than the rule. 

Many authors, beginning with Apostol in the foundational paper [1], use 
transcendental methods when studying the Apostol-Bernoulli polynomials, 
due to their relation with the Lerch transcendent <I>, defined by analytic 
continuation of the series 

$(A, s, a) = 7- — , 

where a 7^ 0,-1,-2,... and either |A| < l,s € C or |A| = l,Res > 1 
guarantees convergence (we use A as a variable in order to maintain a unified 
notation). One has 

(3) Bn{a;\) = -n^{\l-n,a) 

on various domains of analytic continuation, and this relation is often ex- 
ploited to obtain identities satisfied by Bn as special cases of identities satis- 
fied by <I>. The Fourier series of Bn{x; A) is an instance of this (see Section 2). 

While certainly not denying the great merit of this approach, it can be 
overkill, adding unnecessary effort to the study of what is, after all, a poly- 
nomial family. In this paper we wish to bring to light that the algebraic 
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properties of the Apostol-Bernoulli polynomials and basic Fourier analysis 
can get us equally far in describing their asymptotic behavior. We feel this 
is interesting in its own right, more so as the algebraic properties involved 
are precisely those of the so-called "umbral" variety and hence the method 
probably works in a wider context. For instance, this point of view yields 
an overlooked elementary proof of the Fourier expansion (4) itself. 

The same methods used to study the Apostol-Bernoulli polynomials may 
also be applied to the Apostol-Euler polynomials £n{x; A) introduced by Luo 
and Srivastava (see [6, 8]). In Section 8, using a relation which apparently 
has not been previously remarked on in the literature (see Lemma 2), we 
will show that the analogous results for £n Sire consequences of those for Bn', 
hence we concentrate on the latter and summarize the results for the former. 

2. The Fourier series of Bn{x;\) 

In [7, Theorem 2.1, p. 3], it is proved that the Fourier series of Bnix; X) 
for any A € C with A 7^ is 

(4) 8„(.;A) = -M.;A)-3^ Y. p.ft-logA).- 

fc6Z\{0} 

where 6nix; A) = or iog""l according as A = 1 or A 7^ 1. This expansion 
is valid for < x < 1 when n > 2 and for < x < 1 when n = 1. It yields 
the Fourier series of the Bernoulli polynomials, due to Hurwitz (1890), as 
the special case A = 1. 

The expansion (4) is proved from scratch in [7] starting from the gen- 
erating function (1) and applying the Lipschitz summation formula (which 
can be derived in turn from the Poisson summation formula). It may also 
be obtained immediately via the relation (3) to the Lerch transcendent by 
specializing the following series expansion for that can be found in [4, 
formula 1.11 (6), p. 28]: 

(5) $(A, s, x) = A^^T(1 - s) "^{2mk - log Xy~^e^^'^'', 

fcez 

for < x < 1, Res < 0, A ^ (-00, 0]. 

As we have mentioned in the introduction, our approach does not require 
the use of the Lerch transcendent. In this vein, we present an elementary 
proof of (4), using only the algebraic properties of the Apostol-Bernoulli 
polynomials and the real Riemann integral. The main ingredient is the fact 
that {Bn{x]X)} is an Appell sequence for fixed A; namely, it satisfies the 
derivative relation 

(6) B'nix; A) = nBn-iix; A), 

as can be easily deduced from the form of the generating function (1). One 
merely needs to observe that (4) is equivalent to finding the Fourier coeffi- 
cients of X I—)- X^Bn{x; A) for fixed A. 

Proposition 1. For any A € C, X ^ 0,1, k £ Z and n € N, ife have 
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Proof. Fixing k, we use induction on n. One easily checks the case n = 1, 
noting that Bi{x; A) = 1/(A — 1). Let n > 2. Assuming (7) is true for n — 1, 
integrate it by parts, using the derivative formula (6) to obtain 

"1 ^ n! 



/ X^Bnix; X)e-^^'''^dx = AS„(1; A) - e„(0; A) 
Jo 



{27rik - log A)"' 

The proof is concluded by checking that XBn{^', X) = 5„(0;A) for n > 2. 
This follows from (1) noting that Xg{l,X,z) — g{0,X,z) = z (see also [1, 
formula 3.5, p. 165]). □ 

Remark 1. By uniqueness of the Fourier series, (7) could be used to define 
the Apostol-Bernoulli polynomials, instead of using the generating func- 
tion (1). Indeed, it is not hard to work "backwards" to deduce (1) from 
(7) and also to prove the basic algebraic and differential properties satisfied 
by the family directly from the Fourier series, thus providing an alternative 
approach to the theory. This applies also in particular to the Bernoulli and 
Euler polynomials. 

For our purposes, it is useful to rewrite (4) in the alternative form 



(8) ;^^n(^;^) = -E 



ni — ' a" 

e5 



where S is the set of poles of the generating function g{x^ A, z), namely 

if A 7^1, 



(9) S = {ok : k G Z}, Ofc = 27rik — log A, Z 



Z\{0} ifA = l. 



Note that the derivative relation (6) is immediate from (8) by differentiating 
term by term. 

This representation also suggests a proof of (8) from (1) via the calculus 
of residues. Indeed, a~"e"'^ = Kes{z~^'^~^^^g{x,X,z),a) for a G 5. The 
left hand side of (8) is by definition the coefficient of in the power series 
expansion of g around 0, so it is equal to Kes{z~^^~^^^ g{x, A, z),0). Hence (8) 
is equivalent to X^^gc Kes{z~^^~^^^ g{x, A, z), a) = 0, which can be proved by 
showing that the integral <f^^ z~^^~^^^ g{x, X, z) dz tends to over a suitable 
increasing sequence of circles Cn- 

This is essentially the proof in [4, Section 1.11] of the Fourier expansion 
of the Lerch transcendent (5). The recent paper [2] applies this method to 
obtain the Fourier series of the Apostol-Bernoulli, Apostol-Euler, and also 
the Apostol-Genocchi polynomials. Even though it uses no more than basic 
complex analysis, all things considered, after filling in all the details (for 
example, the case n = 1 must be estimated separately from n > 2), the 
proof of Proposition 1 given above is simpler. Nevertheless, the form (8) 
of the Fourier expansion is particularly well-suited for studying asymptotic 
approximations to Bn{z\ A). 

3. Approximations to Bn{X) 

To obtain approximation results from (8) we need to order the set of poles 
S of the generating function (1) by order of magnitude. 

Lemma 1. Let aj^ = 2TTik — log A with /c G Z, A G C, A ^ 0. 
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(a) //ImA > 0, then for k >1, we have 

< |ao| < |ai| < |a„i| < • • • < ja^l < ja_fc| < • • • 

(b) //ImA < 0, then for k >1, we have 

< |ao| < |a_i| < joij < • • • < |a_fc| < \ak\ < • • • 

(c) // A > 0, then for k >1, we have 

\ao\ < = < • • • < \ak\ = |a-fcj < • • • 

(note that qq = if and only if X = 1, in which case it is not a pole, 
hence is excluded in this chain). 

(d) // A < 0, then for k >0, we have 

< jaoj = |aij < |a_i[ = [02! < • • • < |a_fc| = |afc+i| < • • • 
In addition, \ah\ > 27r(|A;| — ^) if \k\ > 1. 

Proof. Let ^ = so that = 2Txi{k — i). Since we use the principal 
branch of the logarithm, this maps A e C \ {0} to the strip -\ <Re^ <\, 
with A > corresponding to Re^ = and A < to Re ^ = ^. Since 
\ak\ = 27r|A;— ,^1, the above chains are readily verified by considering = 
{x — Re^)^ + (Im^)^ for real x. For k G Z, \k\ > 1, we have |A; — ^| > 
\k - Re^l > \k\ - \ and hence \ak\ > 27r(|fc| - \). □ 

By definition, in an asymptotic expansion, the n + 1st term should be 
infinitesimal with respect to the reth term. Thus we have to consider par- 
tial sums of the Fourier series along the chains detailed in Lemma 1 which 
truncate the chain at a link where there is strict inequality. In other words, 
we consider only those finite subsets of poles F C S satisfying 

max{|a| : a € F} < min{|a| : a € S \ F} *== fi. 

Saying that (8) gives an asymptotic expansion means that the remaining 
tail is of the order of 

Let us establish some notation. As in (9), we let Z = Z if A 7^ 1 and 
Z = Z \ {0} if A = 1. Letting = 2'Kik — log A, we define 



Zn 



(10) 



Z„ 



{k£Z 
{k£Z 

{kez 



Fm = {ak ■ k e Zm}, 



\k\ < m}, 

\k\ < m} U {m + 1}, F+ = {ak:ke Z+}, 

|A;| <m}U{-(m+l)}, F~ = {ok : k e Z-} 

for an integer m > 0, noting that for A = 1, Zq is empty. In accordance 
with Lemma 1, we have the following possible choices for F, along with the 
corresponding value(s) of ^ = min{|a| : a G S \ F}: 
(11) 

|a-m+l|, |a-(m+l) I) 
l^"{m+l)l' |Om+l|, 

. J"'-{m+l)l' 

Note that = |a±(m+i)l ^ 27r(m + ^) > vr in each of these cases. 

We begin our study of (8) by proving that it gives an asymptotic expan- 
sion for the Apostol-Bernoulli numbers ;B„(0; A) = -B„(A). 



( TP P+ 

Fm ; Fjy^ , 



l^ 



F+ 



if Im A > 0, 
if Im A < 0, 

if A > (?n > if A = 1) 
if A < 0. 
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Proposition 2. Given A G C, A 7^ 0, let F be a finite subset of the set of 
poles S of the generating function (1) ofBn{x;X) satisfying 

max{|a| : a G F} < mm{|a| : a £ S \ F} = 

For all integers n > 2, we have 

(12) ^ = -Y.^+oi,~n, 

where the constant implicit in the order term depends only on A and F. 

In this sense then, using the appropriate approximating sums over the sets 
-F, the Fourier series (8) of Bn{x; A) at x = is an asymptotic expansion for 
the Apostol-Bernouhi numbers as n 00. 

Proof. Relabel the set of poles in increasing order of magnitude as |ao| < 
|cKi| < ••• < I^mI < •••• The estimate \ak\ > 27r(|A;| — ^) from Lemma 1 
shows that a^*^ is absolutely convergent for n > 2. For any M > 0, we 
have 



00 

-— = - — - x; 

k=M+l l"'^' f,^^j^^ 



1 1 



< 

|aA/+i| 



Cm, A 
|aA/+i| 



^ E 

k=M+l 

where ca/,a is a constant depending only on M and A. This applies to any 
tail, in particular to the tail S \F. □ 

Remark 2. We can say more about the constant in the above proof by 
considering the Hurwitz zeta function for real values of the parameters. For 
example, for F = Fm = {k £ 7^ : \k\ < m} (excluding if A = 1) and 
ImA > 0, for which the next pole in order of magnitude is |am+i|, if we 
denote the constant by Cm,x and let ^ = 3!^, we have, by Lemma 1, 



Cm,A = \am+l P I 12 = |"l + 1 - 



|fc|>m+l ' ' |fc|>m+l ' ' 

oo ^ 

<|m + l-eP ^ -2 =2|m + l-CpC(2,m+i). 

|fc|>m+l (rl ~ 2) 

Simply comparing the sum C(o", q) = YlT=oi^ + for cj > 1, g > with 
the corresponding integral yields 



o- - 1, 

which applied to the above estimate for Cm,\ gives 

m + 



Cm,x < 2|m + 1 - 1 + f -2 • 



Since limm-5.oo("i+l— ^)/(m+^) = 1, this shows that for fixed m and n ^ 0, 
Cm, A is bounded independent of m and we can then replace it with a constant 
cx depending only on A. Thus the tail can be estimated by cx\am+i\~^ for 
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fixed m and n ^ 0. In general, for any n > 2, we still have an estimate 
of the form cx\m + 1 — ^| and hence the tail can always be estimated by 



4. Approximations to I3n{z; A) on the complex plane 

We now come to the central result of this paper, that the Fourier series 
(8) of Bn{x; A) for x G [0, 1] extends to the complex plane as an asymptotic 
expansion for Bn{z;X) given z G C. We state this in a more precise form, 
since we also obtain information on the implicit constants in the asymptotic 
approximation. 

Theorem 3. Given A G C, A 7^ 0, let F be a finite subset of the set of poles 
S of the generating function (1) ofBn{x;X) satisfying 

max{|o| : a G F} < min{ja| : a G S \ F} /i. 
For all integers n > 2, we have, uniformly for z in a compact subset K ofC, 



(13) ^ = -Vi^ + 



where the constant implicit in the order term depends on X,F and K. 

In fact, for n ^ 0, the order constant may be taken equal to the value 
of the constant for the Apostol- Bernoulli numbers, corresponding to z = 0, 
thus eliminating its dependence on K (how large n has to be of course still 
depends on the compact set K). 

Proof. We will use Proposition 2, which is the case 2: = 0, and the "binomial 
formula" 

(14) Bn{x + y;X) = Y,{'t)^n~k[x;\)y'' 

which may be proved directly from the generating function (1) or from the 
derivative relation (6). It is well-known that for a given polynomial family, 
in fact (14) is equivalent to (6) (for fixed x, (14) is the Taylor expansion of 
Bnix + y; A) in powers of y). Both serve to define the notion of an Appell 
sequence. For z G C, writing z = z + 0, (14) and Proposition 2 yield 

Bn{z;X) ^ \ - Bn-kW z^ ^ \ - / ''^ \ n( -n+k\ \ £^ 

n! ^(n-k)\k\ ^\ 2^ a^^^ '] kl 

k=0 ^ ' k=0 \ a£F / 

n ^ u n u 

a6F fc=0 fc=0 

where the implicit constant c is that corresponding to z = and only de- 
pends on F, A. Now, consider the partial sums and tails of the exponential 
series, 

n 00 ^ 

en(w^) = X]-^' e*^{w) = e"" - en{w) = ^ —. 

k=0 ' k=n+l 



L. M. NAVAS, F. J. RUIZ, AND J. L. VARONA 



The first summand above is 



a' 



e*(az) 



To prove the theorem, we need to show that the resulting extra terms have 
the correct order of magnitude. Indeed, 



fc=0 



k=0 



fc=0 



The tails of the exponential series may be estimated using the complex 
version of the Lagrange remainder for Taylor series (which is actually an 
upper bound for the modulus rather than an equality as in the real case): 

L ,|n+l 

K{w)\ < e^"^"^^ + ^^^("'^ = max{Re(u'),0}. 

Since \a\ < fi for all a € -F, we have 



< a e 



\az\_ 



\n+l 



(n + 1)! 



< fie 



I n+1 



(n + 1)! 



so that 



E 



e*Jaz) 



lj.\z\_ 



|n+l 



(n+1) 



(here i^F denotes the number of elements in F) and the latter term is less 
than ce'^l^l/i~" when 



#F- 



< c, 



(n + 1)! 

a condition which certainly holds for n ^ 0, uniformly for z in a compact 
subset K C C. In general, since this sequence is bounded, for any n > 2 we 
obtain a constant independent of n, depending on A, F, K. □ 

Remark 3. Since the point is that we get an asymptotic series in n, the e^'^' 
term could be absorbed into the constant in the order term, but doing this 
hides the numerical behavior of the approximation. Indeed, it is quickly 
apparent in computation that, even for rather small values of j^l, the rapid 
growth of this factor offsets the action of until n is quite large. 

Example 1 (Bernoulli polynomials). As a special case of the theorem, we 
derive Dilcher's results in [3] for the approximation of Bernoulli polynomials. 

Corollary 4. The Bernoulli polynomials satisfy, uniformly on a compact 
subset K of C, the estimates 

(-l)"-i(27r)2'^ 



(15) 



2(2n)! 

(-l)"-^(27r)^"+i 
2(2n + 1)! 



B2n{z) = cos(27rz) + O 



B2n+i{z) =sin(27rz) + 




where the implicit constant depends on K. Moreover, for n ^ this con- 
stant can he made independent of K, equal to the constant for the Bernoulli 
numbers, corresponding to the case z = 0. 
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Proof. The Bernoulli case corresponds to A = 1, where ^ = = and 
the set of poles is S = {Inik : k ^ X,k ^ 0}. Hence Theorem 3, after 
multiplying by (27ri)", states that 

('^^ nl = - ^ "1^+^ {m + ir ■ 

0<|fc|<m \ ^ ^ / 

The result stated above is the case m = 1. □ 



5. Asymptotic behavior of i3„(z;A) 

From now on, assume A 7^ 1. Then the pole set S = {ok = Inik — log A : 
k € Z} contains = — log A and hence the series (13) begins with the term 
_gao2^-n _ (_i)"-i^-2 j^gji j|. jjjg^j^gg ggnse to normalize the series by 

considering the modified expression 
(17) /3„(z;A) = (-i)"-ii^A^0„(z;A). 



nl 



Now (4) becomes 



2iTikx 2Tri 

(18) /3„,(x;A)=V- -— , xG[0,l], A 



^ (1 - Afc)"' ' ' " log A' 

It is straightforward to check that the transformation A = maps the 
region C \ ((-00, 0] U {1}) to the region {A G C : |A - 1| > 1, |A + 1| > 1} 
which is the exterior of the "figure eight" formed by the union of the two 
closed disks |A it 1| < 1, tangent at 0. The negative real axis (— oo,0) is 
mapped to the circumference |A — 1| = 1 minus A = 0. Since |1 — A;A| = 
|1 — k^~^\ = — ^1, where ^ = ^^f^, the ordering of the terms is the 

same as in Lemma 1. 

The asymptotic approximation (13) changes accordingly. 

Proposition 5. Let A G C, A ^ 0, 1. For any integer m >0 we have, for z 

in a compact subset K of C, 

(19) 

I nx 2-Kikz / (2|logA|+27r(m+l))|z|\ 



keim ^ log a; 



_j_ 27r't(m+l) 
log A 



The constant implicit in the order term depends on X,m and K . For n ^ 
it can he taken equal to the constant for the case z = 0, thus making it 
independent of K . Here Im = -^m or are as in (10), with the signs 
chosen according to (11). 

Proof. It is straightforward to check that the general term in the sum changes 
to the expression above, and the order term gets multiplied by lA^log^Aj. 
Estimating |A^| < el^'°S'^l and /i = |a±(m+i)l < 27r(m + 1) + | log A| gives the 
order term. □ 

Remark 4. If we don't incorporate the term A^ into the definition (17) of /3„ 
above, we can drop the "2" from | log A| inside the O term. 
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Corollary 6. For A € C, A 7^ 0, 1 and A ^ (— oo,0), we have, for z in a 
compact subset K of C, 

, W" A / (2|logA|+27r)|^| 

(_l)n-ii^A^E„(z;A) = l + ' 
n! 



mm 



\ -|- 27ri 
log A 



The constant implicit in the order term depends on A and K . However, for 
it can he taken equal to the constant when z = 0. In particular 

lim(-l)"-il^e„(z;A) = A-^ 
uniformly on compact subsets of C 

Proof. This is the case m = above, with the choice of Iq = Zq = {0}. □ 

6. Oscillatory phenomena 

When A G M, A 7^ 0, Lemma 1 shows that there are poles of equal modulus, 
which must be grouped in pairs in the asymptotic approximation. These 
pairs are indexed differently depending on the sign of A, but their behavior 
turns out to be the same. For A > 0, equal moduli poles correspond to pairs 
of integers zizk with A; > 1, and are conjugate: 

a,fc = 2'irik — log A, a_fc = —lirik — log A = o^. 

Writing the poles in polar form, a/; = p^e^'^*"* with = |afc| > and 
Ofc £ [0, 1], this pair contributes to the asymptotic expansion (13) with 

al +al- I' ' +e e 

This simplifies to 

(20) 2A"V^" cos(27r(A;z - na^)). 

Similarly, for A < 0, equal moduli poles correspond to pairs of integers 
{— /c. A; + 1} with A; > 0. It is easy to see that these are also conjugate: 

a-k = —2TTik — log A = — (2A; + 1)111 — log(— A), 

ak+i = 2-ni{k + 1) - log A = {2k + l)Tii - log(-A) 

and hence, now writing aj^+i = pfce^'^*"*^ with pk = |afc+i| > and au G 
[0, 1], the same reasoning shows that this pair contributes to the asymptotic 
expansion with 

(21) 2(-A)"Vfe " cos(7r((2A; + l)z - 2nak)). 

These terms cannot vanish unless z € M, in which case we will write z = x. 
When this is the case, having fixed k and x, the way that the expressions in 
(20) and (21) vary with n depends on the behavior of the sequences kx — nak 
and (2A; + 1)2; — 2nak modulo 1, or equivalently, their fractional parts. 

The study of the fractional parts of /3 — na where a, /? G M are fixed and 
n varies is the problem of real inhomogeneous Diophantine approximation. 
This is a difficult problem with an extensive literature and important ram- 
ifications. Its origins lie in Kronecker's Theorem, which states that if a is 
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irrational, the fractional parts are dense in [0, 1]. In fact, they are equidis- 
tributed modulo 1. Of course if a is rational, the sequence is periodic of 
period the denominator of a. 

As far as we are concerned, then, such oscillatory terms may be considered 
"chaotic" unless is rational, which corresponds to the quotient of the 
conjugate pair of poles afc,afc being a root of unity. 

The simplest appearance of an oscillatory term occurs in the "0th order" 
case of (13) for A G (— oo,0). 

Proposition 7. Let A < 0. Then as n —)■ oo, for z in a compact subset K 
ofC, 

'-^'"'''°^"^-A'g..(.;A) = l+f'°^!t!+"Ve-"' + ( 



n! ' Vl°§ l-^l ~ 



1 -I- 27^^ 
^ log A 



The constant implicit in the order term depends only on A and K . Moreover, 
for n ^ it may be taken equal to the corresponding constant for z = 0, 
thus making it independent of K. 

The part of the approximating terms in parentheses has modulus 1, hence 
is periodic or dense in the unit circle according as its argument is a rational 
or irrational multiple of 2tt. In particular, 

lim f (-')"-' A) - f^iy-VeH = 1 



71— >00 



log |A| — vri 



uniformly on compact subsets of C 



Example 2. It is easy to check that when A < 0, periodic behavior occurs in 
the 0th order approximation if and only if A = —e'^^°^~d for integers d>l 
and k. A special case is A = —1, where we obtain, after simplifying. 



(22) 



2(2n)! 

(_4)n-1^^2n+l 

2(2n + 1)! 



l\n-1^2n 

^2n{z] — 1) = cosvrz + O 



I32n+i{z; -1) = sinvrz + O 




These approximations are essentially those of the Euler polynomials (see 
Section 8). 

Remark 5. An interesting question related to oscillatory behavior is whether 
we can have repeated instances of "good", i.e., periodic, behavior. In other 
words, can we have different pairs of conjugate poles whose quotients are 
roots of unity? When A > 0, considering the Mobius transformation M{z) = 
jz^, which satisfies M(tanz) = e^*^, we have, for k > 1, 



'log A 

and hence the quotient is a root of unity if and only if = tanvr^fc 

with 6 J. € Q satisfying < \9k\ < ^ (this expression cannot be zero). The 
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situation for A < is similar, resulting in 
(2k + 1)tt 

\ I ' = tanvr^fc, Ok G Q, < |0fe| < \. 
log(-A) 

If this happens for two different pairs, corresponding to /c, / > 1, say, then 

tanvr^t. k 2k-\-\ ^ i 

tan TiOi I 21 + \ ^ 

respectively according as A > or A < 0. Without loss of generality one may 
assume both Oi^^Qi > 0. For example, tan ^ / tan ^ = 3 so there is a solution 

pair (A:, 3k) for A = e'^'^^^. However, the ratio k/l = 2 is impossible, as can 
be seen by considering cyclotomic polynomials. What happens in general? 

A similar phenomenon was studied by the authors for Bernoulli polyno- 
mials in [5], namely, when the first conjugate term in the Fourier series fails 
to provide information because it vanishes, one turns to the next, which in 
that case cannot vanish simultaneously with the first. 



7. Successive quotients 



Proposition 8. Let A ^ (0, +oo) and (3n{z; A) as in (17). Then f3n{z; A) 7^ 1 
for n S> and if we denote by e = ±1 the sign of Im A when A ^ R, and 
e = 1 when A € (—00, 0), then we have 



Pn+i{z;\)-l ^ 1 



+ 



2m 
log A 



2-Ki 



n 



00, 



uniformly for z in a compact subset of C In particular 

/3„+i(z;A)-l 1 
nm — — — = TT-^ 

n /3„(2;A)-1 

uniformly on compact subsets of C 

Proof. If A ^ M we may assume without loss of generality that Im A > 0, so 
that e = 1. Applying Proposition 5 with index set / = Zq = {0, 1}, we have 



/3niz;X) = 1 + 



(1 - A)" 



+ 0(|1 + A|-") 



where A = and for simplicity we have subsumed the exponential in |z| 
into the order term. Therefore 



Mz; A) - 1 



1 



^1 - A)" 



^2mz ^ Q 



1 - A 



1 + A 



Since |1 — A| < |1 + A| the e^'^*^ term dominates for n ^ 0, hence I3n{z; A) — 
1 7^ for 77, S> 0. Division of the approximations for n and n + 1 then yields 
the result. 

By Lemma 1 and (11), for A G (—00, 0) we have the same situation, since 
/ = Zq is also a suitable index set. In this case, we have |1 — A| = 1 < 
|1 + A|. □ 
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Remark 6. Compare this result with Corohary 6, which shows that Pniz', A) — 
1 tends to uniformly on compact subsets of C for A ^ (— oo, 0). As we saw 
in Proposition 7, for A < 0, the difference /3„(z;A) — 1 oscillates, yet by 
Proposition 8, its successive quotients approach a limit. 

When A € M"*", we cannot use but rather Zi = {0, 1, —1} as index set. 
This brings into play the oscillatory phenomena mentioned in Section 6, 
which prevent us from obtaining an analogous clean result. However, we 
can still deduce some bounds from (20), or rather a slight modification after 
normalizing. 

Proposition 9. Let /3„ be as in (17). Let A € (0, oo), A 7^ 1, A = 27rz/log A 
and p = |1 — A| = |1 + A|, /x = |1 — 2A|, so that I < p < ^. Let rj = pp^^ , 
so < 7] < 1. Then there is a constant c > independent of n such that 
(23) 

p"-*^! tanh(27rlmz)| - cr?" < 

The estimate holds uniformly in n for z in compact subsets o/C \ M. 

Proof. Let |1 — A| = |1 + A| = p, and p = \1 — 2A|. By Lemma 1, we have 
1 < p < p. Now write 1 - A = pe^'"*" with a £ [0, 1]. Then as in (20), we 
have 

l3n{z; A) - 1 = 2p-" cos 27r(z - na) + 0(p"") 

= 2p-" (cos 2tt{z - na) + 0(77")) . 

In general, if z = x + iy with x, y € M, then | cos = cos^ x + sinh^ y. In 
our case, | cos 27r{z — na) p = cos^ 2Tr{x — na) + sinh^ 27ry. From the remarks 
made in Section 6, unless a is rational, we have no control over the term 
cos 27r(x — na), which will be dense in [—1, 1]. However, if z ^ M then at least 
we have a nontrivial lower bound | cos27r(2: — na)| > | sinh27ry| > 0, which 
eventually dominates the 0((p/i~^)") term. Combining this with the upper 
bound |cos27r(2; — na)| < cosh27ry and forming the successive quotients 
gives the result. □ 

The bounds (23) are sharp. They correspond to the cases when x — na = 
0, |, I modulo 1, where a is defined by w = \\Zf^\ = e^'^*°. When a ^ Q 
the sequence x — na will come arbitrarily close modulo 1 to these values 
infinitely often. 

If we have z = x G M then the cosine term in the approximation Pni-C', A) — 
1 = 2p~"(cos 27r(x — na) + 0{r]'^)) will oscillate densely in [—1, 1] unless a 
is rational. In particular it does not dominate the order term and nothing 
guarantees that /3„(x; A) — 1 is nonzero. Numerical examples where it does 
vanish are easy to find. 

This leaves us with the case z = x £M. where a is rational. Equivalently, 
^ ~ \i~A- \ ~ e^'^*" is a root of unity. Here one can prove the following result. 

Proposition 10. Let j5n be as in (17). Let A G (0,oo), A / 1, A = 
27ri/ log A, p = |1 — A| = |1 + A|, and p = \1 — 2A|, so that 1 < p < p. 
Set rj = pp^^ , so < r] < 1. Suppose that to = p~^(l — A) = e^'^*" is a 
root of unity, i.e. a is rational. Then for z in a compact subset K of C at 



I3n+i{z; A) - 1 
Pn{z;X)-l 



<p ^1 coth(27rImz)| + C7?" 
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positive distance 6 from the exceptional set E{a) C M where cos{x — na) = 
for some n € N, there is a constant c = Ck.x > such that 



(24) p^' ( 



46 



cosh 2iry 



crj 



< 



I3n+i{z; A) - 1 



1 / cosh27r7; 



Proof. Suppose a = ^ with a, d € Z, d > and gcd(a, (i) = 1. Then 
cos27r(x — na) is periodic in n of period (i. It can be checked that E[a) = 
^ + ^Z, ^ + or respectively, according to whether gcd(4a, d) = 1, 2 
or 4. In any case E{a) C ^Z and in fact these three affine lattices are 
disjoint with union ^Z. 

As far as estimation is concerned, we may replace the cosine with the 
function ||^|| denoting the distance from ^ to the nearest integer; namely, we 
have 



2 2a; 



< I cos 27rx| < TT 2x 



In our case, it is easily seen that min„gz ||2x — \ — 2na|| = 2 dist(x, S(q)). 
The rest of the proof is the same as in Proposition 9, using these estimates 
to bound cos 2'k{z — na) below in terms of the distance to the exceptional 
set, rather than using the bound | sinh27ry| as before. □ 



8. The Apostol-Euler Polynomials 

The Apostol-Euler polynomials are a generalization of the Euler polyno- 
mials, introduced by Luo and Srivastava (see [6, 8]). For A G C, A 7^ 1, by 
means of the generating function 

(25) gE{x, A, z) = = V £n{x; A) — , 

Ae^ + 1 ^-^ n\ 

n=0 

which converges for \z\ < \ log(— A)|. For A = 1, we have the classical Euler 
polynomials, f„,(x; 1) = En{x). For A = —1, gE has a simple pole at 2; = 
with residue —2. We can include this case also by defining ^f„(x;A) to 
be, in general, the coefficient of z"" in the Laurent expansion oi qe around 
z = 0, or by redefining qe to be the holomorphic part, i.e. adding —-to 
it when A = — 1. The latter actually makes more sense as far as unifying 
results goes. 

Formula (2.18) of [7] or (37) of [8] gives the relation 

(26) £n{x;X) = ^ (s„+i(x;A) -2"+iS„+i (^;A')) 

between the Apostol-Euler and Apostol-Bernoulli polynomials. However, it 
is easier to use the following relation to transfer results between them. 

Lemma 2. For all A G C, we have 

(27) £n{x; A) = ^0„+i(x; - A). 

n + 1 

In particular, for A = 1, 

n 

(28) Bnix; -1) = --En^iix), n > 1, 
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and for X = —1, 

(29) £n{x; -1) = ^Bn+i{x), n>l. 

n + 1 

Proof. Writing gs for the generating function (1) of the Apostol-Bernoulh 
polynomials, we easily see that gsix, —X, z) = — |g£;(x, A, z), from which 
(27) follows. □ 

Remark 7. The relation (27) shows that the Apostol-Euler and Apostol- 
Bernoulli families are essentially the same. Apparently, (27) is a new ob- 
servation, although we should note that (26) and (27) are related by the 
following "duplication formula" for the Apostol-Bernoulli polynomials: 

(30) rBn (|; A^) = Bnix; A) + ^?„(x; -A). 

Note also that by (28), Example 2 actually describes the asymptotic behavior 
of the Euler polynomials. 

The poles of the generating function gE are the numbers {2k + l)7ri — 
log A for A; E Z. This is also the case for A = — 1, since is a pole of 
gE{x, —1, z). However, the relations (27) and (29) suggest that we exclude 
as a pole. This is consistent with redefining the generating function to be 
the holomorphic part oi gE at 0. Thus we set 



Se 



{{2k + l)Tri-\ogX:k eZ} if A / -1, 
{27rik -.keZ} if A = -1. 



If we write Sb for the pole set of gs, and indicate the dependence on A, this 
means we have the symmetric relation Se{X) = 5'b(— A). 

Now, via (27) we immediately obtain the Fourier expansion of the Apostol- 
Euler polynomials as a special case of that of the Apostol-Bernoulli polyno- 
mials (see also [7]): 

(31) £n{x; A) - ^ 2. ((2fc + l)^i-logA)"+i' 

valid for < x < 1 when n > 1 and for < x < 1 when n = 0. This may be 
rewritten in a form completely analogous to (8) as 

f^^. £n{x;X) _ ^ 

^ ^ 2-n! " ^ a-+i 

which provides further justification for redefining the generating function as 
we have indicated. 

With (32) we now have the same results for the Apostol-Euler polynomials 
as we did for the Apostol-Bernoulli polynomials; indeed, (27) says that they 

are just a special case. We can also think of the poles Ofc *== {2k + l)7ri — 
log A = 27ri — 5 — 'C) with ^ = ^^j-, as "half-integer versions" of the poles 
Ofc = 27ri — log A of gs- In any case. Lemma 1 still describes their ordering, 
and Theorem 3 remains valid with (32) instead of (8). The normalization 
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of £n{z; A) analogous to (17) is given for A 7^ 0, —1 as 
(33) 



Sniz; A) = -/3„+i(z; -A) = (_i)n+i O^gA + ^n)^ e^^'^ X'£n{z; A), 

2 ■ n! 



_ +1 if A > or ImA < 0, 
^ ~ if A < or ImA > 0. 

The sign e is determined by the relation log(— A) = log(A) + eiri. We limit 
ourselves to listing some of the analogous results formulated in terms of 
En rather than Bn, without the specifics about the order of approximation, 
which are of course still valid. For example, the limit in Corollary 6 is now 

(34) lim e„(z; A) = 1, A ^ (0, 00), A / -1, 

n— >oo 

while the oscillating case in Proposition 7 is 

\0gX + Tn\ ^27riz 



n— >oo 



(35) }^Aer.{z;X)-{^^^—-] e^^ ) = 1, A > 0, 



and the limit of quotients in Proposition 8 becomes 

/QftA r £n+i{z;X)-l logA + evri 

36) hm — — = - :, A ^ -00,0 , 

n^oo £n{z; A) — 1 log A — eVTZ 

all holding uniformly on compact subsets of C. 
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